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Abstract. Motivated by recent experiments, which demonstrated lasing and cooling 
of the electromagnetic modes in a resonator coupled to a superconducting qubit, we 
describe the specific mechanisms creating the population inversion, and we study the 
spectral properties of these systems in the lasing state. Different levels of the theoretical 
description, i.e., the semi-classical and the semi-quantum approximation, as well as an 
analysis based on the full Liouville equation are compared. We extend the usual 
quantum optics description to account for strong qubit-resonator coupling and include 
the effects of low-frequency noise. Beyond the lasing transition we find for a singlc- 
or few-qubit system the phase diffusion strength to grow with the coupling strength, 
which in turn deteriorates the lasing state. 
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1. Introduction 

The search for efficient coupling and read-out architectures of scalable solid-state 
quantum computing systems has opened a new field, called "circuit QED" [1]. It is 
the on-chip analogue of quantum optics "cavity QED", with superconducting qubits 
playing the role of (artificial) atoms and an electromagnetic resonator replacing the 
cavity. The resonators can be used to read out the qubit state [2, 3, 4, 5] or to couple 
qubits to perform single- and two-qubit gates [6, 7, 8, 9]. 

Apart from these applications for quantum information processing, circuit QED 
offers the possibility to study effects known from quantum optics in electrical circuits: 
Fock states of the electromagnetic field were created and detected [10], and the 
coherent control of photon propagation via electromagnetically induced transparency 
was shown [11]. In addition, lasing and cooling of the electromagnetic field in 
the resonator has been demonstrated: By creating a population inversion in a 
driven superconducting single-electron transistor (SSET) coupled capacitively to a 
microstripline resonator, Astafiev et al. [12] could excite a lasing peak in the spectrum. 
In another experiment, Grajcar et al. [13] coupled a driven flux qubit to a low-frequency 
LC resonator and observed both cooling and a tendency towards lasing via the so-called 
Sisyphus mechanism. 

In contrast to conventional lasers where many atoms are coupled weakly to the 
light field in a Fabry-Perot cavity, in the micromasers realized, e.g., in Refs. [12, 13] a 
single superconducting qubit is coupled strongly to the microwave field in the resonator. 
Compared to conventional lasers one expects for single-atom lasers a lower intensity of 
the radiation but stronger fluctuation effects. Specifically, the quantum fluctuations of 
the photon number associated with spontaneous emission, which are known to lead to 
the phase diffusion of the laser field [14], have more pronounced consequences. As a 
result, even in the lasing state, phase coherence is lost after a characteristic time r^, 
which sets a limit on the linewidth of the laser radiation and thus on the visibility of the 
lasing signal. Phase diffusion was observed experimentally in a single-qubit maser in 
Ref. [12]; the dependence of the phase diffusion on the coupling strength was analyzed 
theoretically by the present authors in Ref. [15]. 

In the present work we analyze static and spectral properties of single- and few-qubit 
lasers, focusing on the regime of strong qubit-resonator coupling realized in most circuit 
QED experiments. Using a Master equation approach we analyze the consequences of 
qubit-field correlations. In the strong coupling regime they have significant quantitative 
effects on the laser line-width. In addition for the case of few-qubit lasing, we evaluate 
the corrections due to qubit-qubit correlations. They are strongest at the transition 
to the lasing regime but yield only small corrections to the power spectrum. In the 
frame of the so-called "semi- quantum" approximation [16] we describe the qualitative 
differences between multi-atom lasers and superconducting micromasers; specifically we 
analyze the scaling of the lasing transition and diffusion constant with the number of 
atoms. 
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The paper is organized as follows. We start discussing in the following Section the 
two experimental realizations of superconducting micromasers reported in Refs. [12] and 
[2, 13], respectively. In particular we describe how the population inversion in the qubit 
is created in the two examples. In Section 3, we formulate the theoretical model and 
derive the dynamical equations for the micromaser using a master equation approach. 
In Section 4 we review the theory of lasing, paying attention to effects which are usually 
ignored for conventional lasers but are prominent in single- or few-atom lasers in the 
strong coupling regime. We introduce the different approximation schemes. Static 
properties of single-qubit lasers are presented, such as the average photon number and 
the qubit-field and qubit-qubit correlations. We show explicitly that due to spontaneous 
emission in few-qubit lasers the sharp lasing threshold is replaced by a smooth, but still 
well localized transition to the lasing regime. 

Next, in Section 5, we analyze the spectral properties of superconducting 
micromasers. We discuss the effects of correlations between qubit and resonator on the 
phase diffusion process and we show that, beyond the lasing threshold, the linewidth 
grows with increasing coupling strength, thus deteriorating the lasing state. Finally, in 
section 6 we analyze the scaling of the photon number and of the diffusion constant 
with the number of atoms. In addition we demonstrate how low-frequency noise leads 
to inhomogeneous broadening of the lasing peak. 

2. Inversion mechanisms in superconducting micromasers 

2.1. The S SET laser 

The "SSET laser" realized by Astafiev et al. [12] consists of an SSET coupled 
capacitively to a microstripline resonator, as shown in Fig. la. The properties of the 
coupled system and the specific form of the Hamiltonian will be analyzed further in 
later Sections and in Appendix A. In the present Section we describe how a population 
inversion is created in a suitably biased SSET. 

A superconducting single-electron transistor (SSET) consists of two superconduct- 
ing leads coupled by tunnel junctions to a superconducting island. A gate voltage U 
shifts the electrostatic energy of the island and controls, together with the bias voltage 
V, the current through the device. The Josephson coupling, Ej, allowing for coherent 
Cooper pair tunneling through the junctions, is weak compared to the superconducting 
energy gap A sc and to the charging energy of the island, E c = e 2 /2C, C being the 
total island's capacitance. In addition, quasiparticles can tunnel incoherently (with rate 
oc V/eR, where R is the resistance of the tunnel junction) when the energy difference 
between initial and final states is sufficient to create a quasiparticle excitation, i.e., when 
it exceeds twice the gap, \AE\ > 2A SC . We denote by N the number of excess charges 
on the island; it changes by ±1 in a single-electron tunneling process and by ±2 in a 
Cooper pair tunneling event. At low temperatures for the conditions realized experi- 
mentally the number of accessible charge states of the island is strongly reduced. For 
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the further calculations we can restrict our analysis to N = 0, 1, 2. 

We assume the SSET to be tuned close to the Josephson quasiparticle (JQP) cycle, 
where the current is transported by a combination of Cooper pair tunneling through one 
junction and two consecutive quasiparticle tunneling events through the other junction. 
The parameters of the junctions are chosen asymmetrically. By changing the transport 
voltage V and the gate voltage U, we can tune to a situation, where resonant Cooper pair 
tunneling is strong across, say, the lower junction, while quasiparticle tunneling is strong 
across the upper one. Specifically, when the normalized "gate charge" Nq = CqU /e is 
approximately 1, Nq 1, the charge states \N = 0) and \N — 2) are near degeneracy 
with respect to coherent Cooper pair tunneling across the lower junction. Hence the 
eigenstates are 

| |) = cos I |JV = 0) + sin | \N = 2) 

[ ^> = sin I [AT = 0> - cos | [AT = 2) (1) 

where tan£ = Ej/e c h with e c h = 4:(Nq — l)Ec- Quasiparticle tunneling across the upper 
junction leads to transitions between the states | j) and | j) and the odd charge state 




Figure 1. a) The SSET consists of two superconducting leads coupled to an island via 
Josephson junctions with capacitance C\j and Cd- A transport voltage V is applied, 
and the electrostatic energy of the island is tuned by the gate voltage U via the 
capacitance Cq. In addition, the SSET is capacitivcly coupled to an LC-oscillator 
with strength g. 

b) Energies of the charge states \N = 0) and \N = 2) (dashed lines) and the eigenstates 
| t) and | |) (full lines). The energy of the odd charge state \N = 1) may be far from 
the other ones and is drawn at an arbitrary position. If the charge states |iV = 0) and 
\N = 2) are not close to degeneracy, Cooper pair tunneling is suppressed. For Nq > 1 
we have cos | > sin|, and the dominant quasiparticle transitions lead from | J.) to 
| T) ; as indicated by the dashed arrows. The capacitive coupling to the LC-oscillator 
creates an additional coupling between the states | j) and | J.), indicated by the vertical 
arrow. 
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Ji h 

Figure 2. In the setup of Ref. [2] an externally driven three-junction flux qubit is 
coupled inductively to an LC oscillator. 

N = 1). The transition rates are 




The dependence on the relevant matrix elements and the energy gain eV can be lumped 
into the function /(V), which is the normal current through the junction at voltage V. 
Here we can assume that the relevant energy scale for each tunnel event is the applied 
voltage and neglect the smaller change of the energy of the island. 

By choosing £ such that cos | > sin |, we can create a population inversion. In this 
case, the quasiparticle tunneling processes (2) leading from | |) via 1 1) to | f) become 
stronger than the processes in opposite direction (see fig. lb)). From the transition 
rates (2) we readily obtain the bare population inversion in the system, 

D = ^^ = cos£ (3) 
Ptt + Pll 

where pa is the population of the state \i) with i =], {. 
2.2. The superconducting dressed- state laser 

In this subsection we consider systems as investigated in Ref. [2, 13] and shown in 
Fig. 2. Here, a flux qubit is strongly driven by ac-fields to perform Rabi oscillations. It 
is further coupled to a low-frequency LC-oscillator. In the strongly driven situation the 
physics is most conveniently described in the basis of "dressed states" in the rotating 
frame [17]. The transformation to dressed states modifies the relaxation, excitation and 
decoherence rates as compared to the standard results [18, 19]. As a result, for blue 
detuning of the driving frequency compared to the resonant frequency a population 
inversion is produced in the dressed state basis, which in turn can lead to lasing [20]. 

To illustrate these effects we first consider the driven qubit (ignoring the coupling 
to the resonator) coupled to a bath observable X, 

H = - \ AE o z + hQ m cos (u d t) a x 
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In the absence of driving, f2 R0 — 0, and for regular (i.e., smooth as function of the 
frequency) power spectra of the fluctuating bath observables we can proceed using 
Golden rule type arguments [18, 19]. The transverse noise, coupling to a x and a y , 
is responsible for relaxation and excitation processes with rates 
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m (X ),=-ab , (5) 
while longitudinal noise, coupling to a z , produces a pure dephasing with rate 

r; = l ^s x (u = o) . (6) 

Here b± = b x + ib y , and we introduced the ordered correlation function {X 2 ) w = 
f dt e tuJt {X(t)X(0)), as well as the power spectrum, i.e., the symmetrized correlation 
function, S x (ou) = ((X 2 ) w + (X 2 )_ U! )/2. The rates (5) and (6) also define the relaxation 
rate 1/Ti = 11 = 11 + 11 and the total dephasing rate 1/T 2 = Y v = r\/2 + T* which 
appear in the Bloch equations for the qubit. 

To account for the driving with frequency uj d it is convenient to transform to the 
rotating frame via a unitary transformation U r = exp (— iu d a z t/2). Within rotating- 
wave approximation (RWA) the transformed Hamiltonian reduces to 

H = ^MIro cr x + ^h5u a z 

- \ [b z a z + b^a. + b* ± e- i ^ t a + ] X + # bath , (7) 

with detuning 5oo = uo d — AE/h. The RWA cannot be used in the second line of (7) since 
the fluctuations X contain potentially frequencies close to ±u d , which can compensate 
fast oscillations. Diagonalizing the first two terms of (7) one obtains 

H = -hQ'R.O'z + #bath 



^b z + C ^(b* ± e-^ + b ± e^) 



o z X (8) 



rrsin/3 + 1,,, . . sin/3 — 1, • , cos/3, i ~i 

- \ b\ e~ luJdt + b ± e luJdt b z a + X + h.c. [ , 

v - 4 4 2 - J 

where the full Rabi frequency is £7r = a/^ro + Suj 2 , and the detuning determines the 
parameter f3 via 

t&n(3 = 5u/tt R0 . (9) 

From here Golden-rule arguments lead to the relaxation and excitation rates in the 
rotating frame as well as the "pure" dephasing rate [21] 

f^J^cos 2 /3(X 2 ) nR 
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1 - sin/?) 2 (X 2 )^^ + (1 + sin/3) 2 (X 2 ) 



^cos 2 /3(X 2 >„, R 



16ft 2 



(1 - sin/3) 2 (X 2 )_^ R + (1 + sin^) 2 {X\ d _n R 



^sm 2 [3S x (u = 0) + 



Ah 2 



cos 2 S x (uj d ). 



;io) 



We note the effect of the frequency mixing. In addition, due to the diagonalization the 
effects of longitudinal and transverse noise on relaxation and decoherence get mixed. 
We further note that the rates also depend on the fluctuations' power spectrum at the 
Rabi frequency, (X 2 )±q r . 

For a sufficiently regular power spectrum of the fluctuations at frequencies u m 
±AE/h we can ignore the effect of detuning and the small shifts by ±J1r as compared 
to the high frequency uj d as AE/h. We further assume that JIr <C kT/h. In this case 
we find the simple relations 



(1 + sin/3) 5 
4 

(l-sin/3) 2 



r, + 



4 

(1 + sin/3) s 



r T + -cos 2 /3r, , 

r T + \ cos 2 (3Y V , 



(12) 



f;=sin 2 /3r; + ^(r i + r T ), 

where the rates in the lab frame are given by Eqs. (5,6) and the new rate 
T v = ±- 2 b 2 z S x {n K ) 

depends on the power spectrum at the Rabi frequency. 

To proceed we concentrate on the most relevant regime. At low temperatures, 
k B T <C AE as hw d) we can neglect 1^ as it is exponentially small. We also assume that 
can be neglected as compared to Y^, which is justified, e.g., when the qubit is tuned 
close to the symmetry point where b z <C |6j_|. Since the rate T v depends on the noise 
power spectrum at the frequency Qr, which is usually higher that the frequency range 
of the 1/ / noise, the latter does not change the situation. Thus we neglect Y v and we 
are left with 

~ (It sin/3) 2 - cos 2 /? 

r l/T ~ 4 T l » r ^ ~ r l • ( 13 ) 

The ratio of up- and down-transitions depends on the detuning and can be expressed 
by an effective temperature. Right on resonance, where (3 = 0, we have Tf = f^, 
corresponding to infinite temperature or a classical drive. For "blue" detuning, (3 > 0, 
we find If > r_|_, i.e., negative temperature. This leads to a population inversion of the 
qubit, which is the basis for the lasing behavior which will be described below. 

In a more careful analysis, paying attention to the small frequency shifts by ±f2 R , 
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Figure 3. Relaxation rates in the basis of the dressed states. The left staircase 
denotes the eigenstates of the undriven qubit, | [) and the (quantized) driving 
field, |m), before the driving is switched on, i.e., for JIro = 0. The Hamiltonian in 
this basis is obtained from (4) by replacing MIro cos to^t with \(S + d), where A is 
the coupling between the qubit and driving field, and d, dJ are the annihilation and 
creation operators of the driving field. The right staircase stands for the dressed 
states of a driven qubit near resonance, obtained by diagonalizing the corresponding 
2x2 Hamiltonian. This yields \g) = cos - f J \m, f) + sin - £\ \ m + 1, |) 

and |e) = — sin — |j |m, f) + cos (j — |m + 1, J.). The bare Rabi frequency 

is SIrq w where fh is the average photon number in the coherent (classical) 

driving field. In the lab frame at low temperature only the relaxation rate T± needs to 
be considered. However, in the dressed states basis the dominant rate T| leads from 
the state \g) to the state |e), thus creating a population inversion. 



we obtain for (3 = 



(14) 



For example, for Ohmic noise and low bath temperature this reduces to r^/T^ « 
l+2Q^/ud, which corresponds to an effective temperature of order 2hu>d/kB ~ 2AE/k^, 
which by assumption is high but finite. The infinite temperature threshold is crossed 
toward negative temperatures at weak blue detuning when the condition 

(1 + -/»'.. 1 + 2k (15) 



(1-sin/?) 2 u d 

is satisfied. We note that all qualitative features are well reproduced by the 
approximation (13). 
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To illustrate the calculations outlined above and the mechanism creating the 
population inversion for blue detuning we show in Fig. 3 the level structure, i.e., the 
formation of dressed states, of a near- resonantly driven qubit. For the purpose of the 
present discussion we assume that also the driving field is quantized. This level structure 
was described first by Mollow [17]. The picture also illustrates how for blue detuning a 
pure relaxation process, T^, in the laboratory frame predominantly leads to an excitation 
process, If, in the rotating frame, thus creating a population inversion in the basis of 
"dressed states" . 

If this effective inverted two-state system is coupled to an oscillator, a lasing state 
is induced. In Ref. [20] it was proposed to couple the oscillator to the dressed states 
belonging to the neighboring doublets (see Fig. 3). Then, to be in resonance with the 
pair of dressed states with population inversion the oscillator frequency should satisfy 
ujq = uJd + Qn- As ood ~ AE and Q-r, *C AE this can work for a high-frequency resonator 
approximately in resonance with the qubit, c^o ~ AE. In contrast, in Ref. [22] a different 
situation was considered where the oscillator was coupled to the dressed states belonging 
to the same doublet. The resonance condition then reads ojq = JIr, and the lasing can be 
reached for an oscillator much slower than the qubit, ojq <C AE, which is the situation 
realized in Ref. [2]. An additional complication arises at the symmetry point of the 
qubit, since there the single-photon coupling between the oscillator and the doublet 
of the dressed states vanishes. Then, two-photon processes become relevant with the 
resonance condition 2u; = [22]. 

3. Modelling the single- or few-qubit laser 

We consider a single- mode quantum resonator coupled to N a qubits (labelled by /i). 
In the absence of dissipation, in the rotating wave approximation, the dynamics of the 
system is described by the Tavis-Cummings Hamiltonian [23]: 



Here we introduced, apart from the photon annihilation and creation operators, a and 
a\ the Pauli matrices acting on the single-qubit eigenstates = | j^) (|^| — (J, M |, 
° + = It//) UmI> an d a - = II//) ( Tm I • Including both resonator and qubit dissipation the 
total Hamiltonian becomes 



Dissipation is modeled by assuming that the oscillator and the qubits interact with noise 
operators, X a and X£, X+, X^, belonging to independent baths with Hamiltonian ff bath 
in thermal equilibrium [24]. The noise coupling longitudinally to the qubits, X£a%, is 
responsible for the qubits' pure dephasing. 

In this section and beyond we do not describe anymore the detailed mechanism 
creating the population inversion in the qubits, which is necessary to obtain lasing. 




(16) 



H = H TC + (a + a^)X a + ^ (X£a% + X%a% + Xtat) + if, 



bath- 



(17) 
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Rather we introduce it by assuming that the effective temperature fixing the ratio of 
excitation and relaxation rates of the qubits is negative. (In the same spirit the transition 
rates V appearing below are those of the effective two-level system, even if they refer 
to transition between dressed states, for which the rates were denoted above by V.) 
Possible deviations from the Tavis-Cummings oscillator-qubit coupling used in Eq. (17) 
are discussed in Appendix A. 

The dynamics of a single- or few-qubit laser can be analyzed in the frame of a 
master equation approach, as discussed by several authors [22, 26, 27, 28, 29]. In the 
Schrodinger picture the master equation for the reduced density matrix p of the qubits 
and the oscillator reads 

p = -^[H TC ,p]+L Q p + L R p. (18) 

The Liouville operators L R and Lq describe the resonator's and qubits' dissipative 
processes. For Markovian processes it is sufficient to approximate them by Lindblad 
forms, 

Z/r p = — [(iV t h + 1) (2apa) — (Jap — pa) a) 

+ N th (2a f pa - aaV - paa ] )] (19) 



and 



L QP= 



r 



(20) 



The dissipative evolution of the system depends on the excitation, relaxation, and pure 
dephasing rates of the qubits, , 1^ and T* , as well as on the bare damping rate of the 
resonator, k, and on the thermal photon number N^. 

For later purposes, we also introduce the rate Ti = + Tf, which is the sum of 
excitation and relaxation rates, and = Ti/2 + T* , the total dephasing rate, incl. the 
"pure dephasing" due to longitudinal noise described by T* . In contrast to relaxation 
and excitation processes, pure dephasing arises due to processes with no energy exchange 
between qubit and environment and thus does not affect the populations of the two qubit 
states. The parameter D = (r-f — T^/Ti denotes the stationary qubit polarization in the 
absence of the resonator. In the present case, since we assume a negative temperature 
of the qubit baths and a population inversion, we have Dq > 0. 

The master equation (18) allows us to determine completely the quantum state of 
the system. However, its full solution is numerically demanding in the experimental 
regime of parameters due to the high number of photons in the resonator (of the order 
of 10 2 or higher for a single-qubit laser). For this reason, we will use, whenever possible, 
different approximation schemes to calculate the physically relevant quantities. 
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4. Approximations and static properties 



To describe the single- or few-qubit laser in the strong coupling regime we start from the 
master equation (18) for the density matrix. In some cases we find that approximate 
analytical results, which are presented in this section, are sufficient. In general, however, 
we rely on a numerical solution. 

From Eq. (18) we obtain the following equations for the average photon number 
(n), the qubit polarization (a%) and the product {o+a), 

j t (a») = -2ig («a> - (atrf)) - I\(K> - D ) , 
j t (n) =igJ2 (««> " " « «n> - N th ) , 



d_ 

dt 



;<a> = (<A - 7 ) «a> - igtfn) - ig 



(21) 



Here we introduced the detuning A = e — uj and the total dephasing rate 7 = + |. 

In the stationary limit, after isolating the correlations between different qubits by 
writing (a+crU) = 5^(1+ + (1 — ^)(c+c!l), we derive from the previous equations 
the following two exact relations between four quantities: the average qubit polarization 
(S z (t)) with S z = jj-J^n^zi ^e photon number (n(t)), and the correlators (nS z ) and 
C QQ = i:E M ^« (7 -)> 



(n> 



N th + 



2g 2 N a 7 



K 



(S z ) = D 



7 2 + A 2 

7 



(S z n) + -((S Z ) + 1) + C QQ 



(S,n) + -((S,) + 1) + C< 



(22) 



ri 7 2 + a 2 

If two of them are known, e.g., from a numerical solution of the master equation, the 
other two can be determined. 



4-1. Semi- quantum model 

Factorizing the correlator, (S z n) ~ (S z )(n), on the right-hand side of Eqs. (22) and 
neglecting the qubit-qubit correlations, Cqq ~ 0, we reproduce results known in 
quantum optics as "semi-quantum model" [16]. This approximation yields a quadratic 
equation for the scaled average photon number (per qubit) n = (n)/N a , 



N n 2iV 2 4k N n 



0, (23) 



which depends on the parameter n = + ■ This equation has always one 

positive solution n > 0. 
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4-2. Semiclassical approach 

Before continuing the analysis of the properties of the semi-quantum solution, for sake 
of comparison, we recall the standard semiclassical results. In this approximation the 
operator a is treated as a classical stochastic variable, a. After adiabatic elimination of 
the qubits' degrees of freedom, i.e., assuming T v k/2, one obtains a classical Langevin 
equation for a, 

- 2 V 1 

a + £(t)- (24) 



a = — 



k _ g 2 N a 
2 T M -iA z 



Here is a classical Langevin force due to thermal noise, (£(£)£*(£')) = K ^th^{t — t'), 
and sf = Dq/(1 + \a\ 2 / (n N a )) denotes the stationary qubits' polarization. 

In order to obtain an expression for the average photon number (n) = (\a\ 2 ) we 
rewrite the Langevin equation as a = — f(\a\ 2 )a+£,(t) and approximate (/(|a| 2 ) • \a\ 2 ) ~ 
f({n)) ■ (n). Thus we arrive at the equation j^{n) = — 2Re{/((n))} • (n) + nN th , from 
which we obtain in the steady state a quadratic equation for the scaled photon number 
n = (n)/N a , 

_ 2 / T 1 D iV th \ N th h 

n + { n °-^-^) n - — = (25) 
In the low-temperature limit, iVth ~ 0, the semiclassical results can be rewritten in the 
simple form, n 2 + (n — TiD /2k) n w 0, with n = (l + fr)- This gives the 

well-known threshold condition D > kY^/ (2g 2 N a ) for the lasing state. 



4-3. Comparison of the different approaches 

Different from the semiclassical picture, the semi-quantum model includes the effects 
of spontaneous emission processes, described by the term proportional to ((S z ) + 1) in 
Eqs. (22). Spontaneous emission is responsible for the linewidth of the lasers, and, as 
noticed in Ref. [27], due to the low photon number, spontaneous emission is especially 
relevant for the dynamics of single-atom lasers. 

To illustrate the effect of spontaneous emission on the lasing transition and at the 
same time the quality of the semi-quantum approximation, we plot the photon number 
as a function of the coupling strength g for N a = 1 (Fig. 4, left panel) and N a — 2 (Fig. 
4, rigth panel). The plots show the semi-quantum, semiclassical and Master equation 
results. We note that the semi-quantum approximation gives results in very good 
agreement with the Master equation. Moreover, both the semi-quantum and Master 
equation solution show a smooth crossover between the normal and the lasing regimes, 
an effect which is due to spontaneous emission. While we cannot define a sharp threshold 
condition, we can still identify, even for a single-atom laser, a well localized transition 
region centered at the threshold coupling predicted by the semiclassical approximation. 

In the left panel of Fig. 4 we also plot the qubit-oscillator correlator, (S z n), and 
the factorized approximation, (S z )(n). For strong coupling, they differ significantly. 
However, as the good agreement between the semi-quantum approximation and the 
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Figure 4. Left panel: average photon number (n) in the resonator and qubit-field 
correlations for a singlc-qubit-laser. The photon number is calculated using the Master 
equation (ME, solid black line), the semi-quantum (SQ, solid orange (light grey) line), 
and the semiclassical approximation (SC, dotted line). The dot-dashed and dashed 
lines show the average values (a z n) and (a z )(n). Right panel: average photon number 
(n) in the resonator and qubit-qubit correlations Cqq (dashed line) for a two-qubit- 
laser. The photon number was calculated using the Master equation (solid black line), 
the semi-quantum (solid orange (light grey) line), and the semiclassical approximation 
(dotted line). The parameters are e = uj , Ti/w = 0.016, r* /uj = 0.004, D = 0.975, 
k/lo = 3 • 10~ 4 , and iV th = 0. 

numerical solution of the Master equation demonstrates, the qubit-field correlations 
have only a weak effect on the average photon number. On the other hand, as we will 
see in the following section, the qubit-field correlations have an important effect on the 
spectral properties of the single-qubit laser. 

The right panel of Fig. 4, shows the qubit-qubit correlations Cqq for a two-qubit- 
laser. Similar as the qubit-field correlations, they are neglected in the semi-quantum 
approximation. Both correlations are maximum at the lasing transition, but decay away 
from this point. The reason is that qubit-field and qubit-qubit correlations scale as 
g 2 /T^, thus they are small for weak coupling. On the other hand, they are proportional 
to the qubit inversion (S z ) and hence vanish rapidly above the transition. 

5. Spectral properties 

In this section, we will study the spectral properties of single-qubit lasers. The 
emission spectrum 0(u) is given by the Fourier transform of the correlation function 
Q(r) = lim t _ 00 (at(t + r)a(t)). As we will see, for typical circuit QED parameters, i.e., 
for strong coupling g, the semi-quantum approximation, in spite of giving a sufficient 
estimate of the stationary photon number, cannot be used for a quantitative study of 
spectral functions. We evaluate the correlation function by performing a time-dependent 
simulation of the master equation (18) using the method described in Ref. [28]. This 



Lasing in Circuit QED 



14 



method is numerically demanding, especially when we consider lasing with more than 
one qubit, N a > 1. We will show that in resonance, e = Uq, the semi-quantum theory 
catches the most qualitative features, both below and above the transition to the lasing 
regime. We will use this method later in Section 6.1 to investigate the scaling of the 
spectral properties with the number of qubits N a . 

5.1. Spectral properties in the semi- quantum theory 

Similarly to Eqs. (21) for the average values, we can derive equations for the laser and 
cross correlation functions, O(r) and G{t) = lim i _ >00 (o" + (t + r)a(t)). Assuming that 
the oscillator damping is much weaker than the qubits' dephasing, k/2 <C T v , which 
is usually satisfied in single-qubit lasing experiments, we obtain a single equation for 
the oscillator correlation function: -^{o){t + r)a(t)) = (iuj — k/2) (a\t + r)a(t)) + 
g 2 /T v (a z )(a^(t + r)a(t)). Thus the semi-quantum theory predicts an exponential decay 
of the correlation function O(r), which corresponds to a Lorentzian shape of the emission 
spectrum, 



5.2. Numerical investigation of the spectral properties 

For the following discussion we focus on the case of a single qubit, N a — 1. In the left 
panel of Fig. 5 we plot the diffusion constant K,d, as function of the coupling strength 
g, covering the whole range from below to above the transition, and compare it to the 
diffusion constant /t|f c , obtained from the semi-quantum theory, i.e., by neglecting the 
qubit-field-correlations (cr z n) — (a z )(n). 

Upon approaching the broadened lasing threshold from the weak coupling side 
we observe the linewidth narrowing characteristic for the lasing transition. However, 
above the transition, the linewidth increases again with growing coupling strength, thus 
deteriorating the lasing state. By comparing the semi-quantum approximation with the 
full solution of the master equation, we observe that qubit-oscillator correlations have 
a significant quantitative effect on the phase diffusion, leading to a reduction of the 
linewidth by roughly a factor 1/2, but they do not change the qualitative conclusions. 

Also in Fig. 5, we note that in the transition region there is an "optimal" value of the 
qubit-oscillator coupling where the height of the spectral line, which is given by the ratio 
(n) I Kd of the photon number and the linewidth, is maximum. This interesting feature 
is due to the fact that in single- and few-qubit lasers far above the lasing transition a 
increase of the coupling has little effect on the saturated photon number, but leads to 
an increase of the incoherent photon emission rate and the linewidth. 




2n d (n) 



(26) 




(27) 
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Figure 5. Left panel: Phase diffusion constant and average photon number (solid 
line) as a function of the coupling strength g for a single-qubit laser. The phase 
diffusion constant is calculated using the numerical simulation (dashed line) and the 
semi-quantum approximation (dotted line). Right panel: Phase diffusion constant 
(dashed line), frequency shift (dotted line) and photon number (solid line) as a function 
of the detuning for a fixed coupling strength g/uj = 0.005. In both plots, we used 
k/uiq = 5 • 10~ 4 for the bare damping rate of the resonator, other parameters as in 
Figure 4. 

When the qubit and the resonator are not in resonance, A ^ 0, the emission 
spectrum is shifted with respect to the natural frequency ujq of the resonator. This 
is shown in the right panel of figure 5 where we plot the average photon number, the 
linewidth and the frequency shift 5uq as functions of the detuning A in the strong 
coupling regime. 

The numerical results for the linewidth shown in this plot differ qualitatively from 
the results presented in our previous paper [15], where we use a factorization scheme 
to obtain analytical expressions for the linewidth. Specifically, moving away from the 
resonance the linewidth is increasing, while the factorization predicted a decrease. On 
the other hand, the approximation based on the factorization yields results very similar 
to the numerical ones right on resonance, as well as in the far off-resonant situation and 
the weak-coupling regime. 

6. Discussion 

6.1. Scaling in the semi-quantum approximation 

As discussed in a quantum optics context in Refs. [27, 28], various properties of 
single qubit masers are due to the fact that in these systems only one artificial atom 
(a microscopic system from a thermodynamical point of view) interacts with the 
electromagnetic radiation. To clarify the main differences between single qubit masers 
and conventional (many atom) lasers, we use the semi-quantum approximation to study 
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the scaling of the average photon number and of the phase diffusion with the number 
of atoms. 

In Fig. 6, we plot the scaled photon number (n)/N a in the transition region and 
K d versus the scaled coupling, g^/Nl, for different values of the number of qubits N a . 
Plotted in these scaled forms, all curves have the same asymptotic behavior, and the 
transition occurs at the same position. As expected, in Fig. 6 (left panel) we observe 
that for low values of N a , there is a smoothening of the lasing transition which is due 
to spontaneous emission processes and disappears in the large N a limit. In Fig. 6 (right 
panel) we show the scaling of the phase diffusion constant. Here the qubits' relaxation 
processes are responsible for the increase of the phase diffusion rate in the case of strong 
qubit-oscillator coupling for small N a . 

6.2. Effect of the low-frequency noise 

The linewidth of order of 0.3MHz observed in Ref. [12] is about one order of magnitude 
larger than what follows from our results (of the order of the Schawlow-Townes 
linewidth). Moreover, in the experiment, the emission spectrum shows a Gaussian 
rather than a Lorentzian shape. Both discrepancies can be explained if we note that the 
qubits' dephasing is mostly due to low-frequency charge noise, which cannot be treated 
within the Markov approximation used in the present analysis. 

However, low-frequency (quasi-static) noise can be taken into account by averaging 
the Lorentzian spectral line over different values of the energy splitting e of the qubit 
[30] , or equivalently, over different values of the detuning A between qubit and oscillator. 
Assuming that these fluctuations are Gaussian distributed, with mean A and width a, 
such that Ti > a ^> Kj, we can neglect in the saturated limit the dependence of and 
(n) on A and assume that the frequency shift 5u depends linearly on the detuning A. 

In the strong coupling regime, the shift of the emission spectrum is given, in a 
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good approximation, by 5u ~ A/c/(2r ¥ ,), which leads to a Gaussian line of width 
o — uk/ (2r v ), where we remark that is the total Markovian dephasing rate. In this 
way, the linewidth observed in the experiment can be reproduced by a reasonable choice 
of a of order of 300 MHz. In the case in which a is larger than r\, the previous formula 
overestimates the linewidth since it does not take into account the decay of (n) below 
the lasing transition. In this case we can still perform the averaging numerically. In 
either case we note that in the presence of low-frequency noise, the linewidth is governed 
not by Kd, but by 5uq. 

7. Conclusions 

We analyzed in detail the static and spectral properties of single- and few-qubit lasers. 
Our main conclusions are: 

- As compared to a conventional laser setup with many atoms, which has a sharp 
transition to the lasing state at a threshold value of the coupling strength (or inversion), 
we find for a single- or few-qubit laser a smeared, but still well defined transition. 
Similarly, the decrease of the phase diffusion strength when approaching the transition, 
i.e., the characteristic linewidth narrowing, is less sharp but still pronounced. 

- Above the lasing trasition we observe for a single- or few-qubit laser a pronounced 
increase of the phase diffusion strength, which leads to a deterioration of the lasing 
state and a reduction of the hight of the laser spectrum. 

- Low-frequency noise strongly affects the linewidth of the lasing peak, leading to 
an inhomogeneous broadening. In comparison, the natural laser linewidth due to 
spontaneous emission is negligible. 

Acknowledgments 

We acknowledge fruitful discussions with O. Astafiev, J. Cole, A. Fedorov and F. 
Hekking. The work is part of the EU 1ST Project EuroSQIP. 

Appendix A. Comment on two-photon processes 

Here we briefly discuss the validity of the Jaynes-Cummings model introduced in Section 
3, when applied to describe the SSET laser of Astafiev et al. [12]. As discussed in 
Section 2.1, the SSET laser, schematically depicted in Fig. 1, consists of a biased 
superconducting island coupled capacitively to a single-mode electrical resonator. Under 
appropriate conditions only two charge states, corresponding to N = 0, 2 are relevant 
to the dynamics of the device. In this basis the hamiltonian of the oscillator and the 
qubit can be written as 

H = \ (e c h t z + Ejt x ) + hw a^a - hg r z (a + a f ) . (A.l) 
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where the operators r x and r z are defined as: r z = (\N = 2) (N = 2| — \N = 0) (N = 0|) 
and t x = (\N = 2) (N = 0| + |JV = 0) (N = 2|). Rotating to the qubit's eigenbasis 
{| |), | I)}, defined by Eqs. (1), we can recast the hamiltonian as follows: 

H — -e<7 z + hu a^a — hg (cos£<t 2 — sin t;cr x ) (a + a^) . (A. 2) 

The angle £ is defined as in Section 2.1, tan£ = and the qubit energy splitting, 
e, depends on the charging and Josephson energies, e c h and Ey. e = a/ + j • ^ n 
order to identify the one- and two-photon coupling strength, we now apply a Schrieffer- 
Wolff transformation U = e lS with 5 = i g0 ^° s ^ o- 2 (a — a' 1 ") and perform a perturbation 
expansion in the parameter go/u} . The transformed Hamiltonian, H = WHU, thus 
becomes 

H ~ -ea z + hjjQO^a + hg\O x (a + a^) + hg^ioy (a 2 — (a^) 2 j • (A. 3) 

Here we neglected terms of order (g /u)o) 3 and introduced the two coupling constants 
gi = —g sin £ and g 2 = sin £ cos £ for one-photon and two-photon transitions, 
respectively. For the parameters used in the experiment the coupling g 2 is roughly two 
orders of magnitude smaller than the one-photon coupling and below the semiclassical 
threshold for the two-photon lasing, g^ r = A//t 2 r v ,/(r 1 Dg) [31]. In the parameter 
regime explored in the experiments, the Hamiltonian used in Eq. (17) gives thus a good 
description of the dynamics of the system. 
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